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SUMMARY 

Step-by-step procedures are described which can be used for the prelimi- 
nary design of fiber composite box beams subjected to combined loadings. These 
procedures include a collection of approximate closed-form equations so that 
all the required calculations can be performed using pocket calculators. 
Included is an illustrative example of a tapered cantilever box beam subjected 
to combined loads. The box beam is designed to satisfy strength, displacement, 
buckling, and frequency requirements. 


INTRODUCTION 

The design of fiber composite structural components requires analysis 
methods and procedures which relate the structural response of the component to 
the specified loading and environmental conditions. Subsequently, the struc- 
tural response is compared to given design criteria for strength, displacement, 
buckling, vibration frequencies, etc. in order to ascertain that the component 
will meet all the design requirements and will perform satisfactorily. 

An important class of structural components that can readily be made using 
fiber composites are box beams. Box beams are generally used to span long dis- 
tances and to resist combined loads. Box beams are the main structural compo- 
nents in aircraft wings. They are made using thin flat/curved laminates, are 
designed to resist the loads primarily through membrane action and are designed 
to have constant or tapered cross sections. In addition, the laminate thick- 
ness for the covers and sides can be different and varied along the span. In 
a previous paper (ref. 1) step-by-step procedures were described for the pre- 
liminary design of composite panels subjected to combined loadings. These pro- 
cedures have since been extended for the preliminary design of composite box 
beams. The objective of this paper is to describe these extended procedures. 

These procedures include a collection of simple equations to expedite the 
various calculations performed during the preliminary design phase. These pro- 
cedures are demonstrated by applying them to a preliminary design of a tapered 
cantilever box beam. The box beam is subjected to combined loads at the free 
end. It is designed to meet strength, displacement, buckling, and frequency 
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requirements. The various steps in these procedures are described in detail 
with ample explanatory notes so that they can be used to aid in the preliminary 
design of built-up composite structural components in general. 

SAMPLE DESIGN 

It is necessary to have as complete a definition of the specific design 
as is possible in order to initiate the preliminary design phase. For the 
illustrative example described herein, this definition consists of the 
following. 

1. Structural Component: 

Cantilever, 3-bay box beam (schematics figs. 1 and 2). 

2. Specified Loads: 

Free end static loads (fig. 1). 

6600 lb vertical; 3300 lb lateral; 100 000 lb in twist moment. 

3. Displacement Limits: 

Tip displacements less than 1.5-percent of length; angle of twist less 
than 1°. 

4. Frequencies: 

Flap greater than 100 cycle/sec, edge greater than 150 cycle/sec; twist 
greater than 450 cycle/sec. 

Local panel frequencies to be greater than box beam global frequencies. 

5. Safety Factor: 

2.0 times specified load. 

6. Composite System: 

As graphite fiber in epoxy matrix (AS/E) about 0.6 fiber volume ratio. 

7. Design Procedure/Requirements: 

Box beam not to exceed displacement limits. 

Laminates in various bays not to exceed ply fiber-controlled strengths at 
design loads or ply matrix controlled strengths at specified loads. 

Composite panels in each bay not to exceed combined stress buckling. 

8. General philosophy on preliminary design of composite box beams: 

Size covers for only the vertical load and add plies for the combined 
loads (lateral and twist moment). 
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Size side walls for only the lateral load and add plies for the combined 
loads (vertical and twist moment). 


STEP-BY-STEP DESIGN PROCEDURE 

Once the design is defined to the extent just outlined, we are' ready to 
design the composite laminates for the covers and the walls of the box beam by 
following the step-by-step design procedure. 


Step 1: Identify Design Variables 

Number of plies, ply orientation and stacking sequence for the composite 
covers and side walls for the three different bays. 


Step 2: Establish Design Loads 

Safety factor times specified loads (fig. 1): 

N cxx = 2 x vertical load (6600 lb) = 13 200 lb 

N C yy = 2 x lateral load (3300 lb) = 6600 lb 

M cxx = 2 x twist moment (100 000 lb in.) = 200 000 lb-in. 


Step 3: 

Obtain composite material properties (ply and ±9 angleply) for AS/E 
from table I and figures 3 and 4. 


Step 4: 

Select laminate configurations for box beam covers and side walls in each 
of the three bays. 

(a) Calculate in-plane membrane loads at the bulkhead locations (figs. 1 
and 2): These loads are calculated by dividing the moment at that section by 

the respective depth and width. 

At span station 0: 

covers: 

N cxx = Pz x 2,/(h x w) 

= 13 200 lb x 60 i n . / ( 1 0 in. x 20 in.) 

= 3960 lb/in. 
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wal Is: 


N(-yy = Py X 2,/(h X W) 

= 6600 lb x 60 in./ 20 in. x 10 in. 
= 1980 lb/in. 


covers : 

^cxy = ±M x /[(w x h) covers + (w x W a 1 1 s ^ 

+ Py /w covers 

= ±200 000 lb-in. / [ < 20 in. x 10 in.) 

+ (10 in. x 20 in. )] 

+ 6600 1 b / < 2 x 20 in.) 

= -335 lb/in. (top cover) 

665 Ib/in. (bottom cover) 

walls: 

^czx = ±M x /[(w x h) covers + (w x wa n$l 

+ Py/h s ide walls 

= ±200 000 lb-in. / C ( 20 in. x 10 in.) 

+ (10 in. x 20 in. )] 

+ 13 200 1 b/(2 x 10 in.) 

= 160 lb/in. (front wall) 

= 1160 lb/in. (back wal 1 ) 

Repeating the calculations for the other span stations and summarizing in 
table form, we have the results shown in table 2. 

(b) Design bottom cover (pressure surface, figs. 1 and 2). This surface 
is in tension. We need to use the longitudinal tensile strength. Number of 
0° plies Nq,o = Design load (N cxx )/(longitudinal tensile strength 
(Sant = 220 000 PSi> x P!y thickness (t& = 0.005 in.)). 

ki _ N CXX +3960 lb/in. , c 

a° = s am t a ■ 220 000 lb/sq in. x 0.005 in. = 3<6 ~ 4< 

Number of ±45° plies N a±4 5 = Design load (N cxy ) x one-half the ratio of the 
ply longitudinal modulus (E^i = 18.5 mpsi) to ±45° composite shear modulus 
(Gqi 2 = 5.8 mpsi ) / 1 ong i tudi nal compressive strength (Soiic = 180 000 psi) x 
ply thickness (t^ = 0.005 in.) 
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N = N cxy x (1/2)(E an /G 912 ) _ 665 1 b/i n (1/2(18.5/5.8) _ . _ 

2±45 s 2,11c X * 2 , 180 000 lb/sq in. x 0.005 in. ' b ~ 

The number of plies is rounded up since plies are available in fixed thick- 
nesses. Number of 90° plies N^gg: Note there is no N C yy design load 

and therefore no 90° plies are needed. However, we will use two plies for 
laminate integrity and improved buckling resistance. Therefore, 


N 290 ~ 2 

Thus, the laminate configuration for the bottom cover is eight plies as follows 
4 at 0 °, 2 at ±45° and 2 at 90°. Note that this is not a symmetric laminate. 

(c) Design top cover (suction surface figs. 1 and 2). This surface is in 
compression. We need to use the longitudinal compression strength to determine 
the number of zero plies. 


"cxx 

^nA 


-3960 lb/in. 

180 000 lb/sq in. x 0.005 in. 


4.4 as 4 


The number of ±45° plies is now determined from the ratio: 
N„ 


N CX y< T 0 P ) 

'2*45 = N cxy (Bottom) x 1,16 plies 


335 

= 505 x 1.16 = 0.58 as 1 -> use two plies for a balanced lamimate 


N^90 = 2 (same reason as for bottom cover) 


Therefore, the top cover is an eight-ply laminate same as bottom cover. 

(d) Design side wall (leading edge (front), figures 1 and 2, and loads 
from table step 4. a) 


cxx 


1980 Ib/in. 


20 


S 211 t x ^2 220 000 11:) /sq x °-°°5 in. 


= 1 .8 * 2 


2±45 


N c,v <1/2>< WW 
S 211c X *2 


160 lb/in.(1/2)(18.5/5.8) 
180 000 lb/in . 2 x 0.005 in. 


= 0.3 plies -» use two plies 


N 290 = 2 (same reason as for covers) 


(e) Design side wall (trailing edge (back) figures 1 and 2 and loads from 
table step 4. a) 
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ki _ N cxx -1980 Ib/in. , , . 

2.0 S^i 1 c x 180 000 lb/sq in. x 0.005 in. ~ L 

N _ N cxy 1160 Ib/in. . 

&±45 S anc x t a ' 180 000 lb/sq in x 0.005 in. = ~ * 

N^go = 2 (same reasons as for covers) 

Therefore the laminate configuration for the sidewall is six plies as follows: 

2 at 0°, 2 at ±45°, and 2 at 90°. Note that this also is not a symmetric 
laminate with respect to bending. 

(f) Select plies at the other span stations. Examining the loads in table 
step 4. a, we see that: (1) N cxx is smaller at all the other span stations, 
and (2) N CX y is maximum at the 60-in. span station. At this point we can 
either calculate the number of plies at each of the span sections or calculate 
the number of plies we need for the maximum N CX y at the 60 in. station and 
use a uniform laminate through the box beam. Using a uniform laminate simpli- 
fies the fabrication procedure but it will increase the weight. Assuming that 
the weight is not critical, we will calculate the number of ±45° plies we need 
at maximum N CX y and use this number throughout the box beam. 

The number of ±45 plies at maximum N CX y (60-in. station) is 

N 0+4 , N cxy E Sl11 /2G ei2 _ 3320 lb in. (1/2X18.5/5.8) - Q 

S aiic t a ~ 180 000 lb/s( i in - x °- 005 1n - " ~ b 

Therefore, the laminate configuration for the box beam is 10 plies as follows: 

4 at ±45°, 4 at 0°, and 2 at 90°. Using the conventional notation, this is 
expressed : 


[±45/0 2 /90] s 

(g) Select the minimum number of plies at the various span stations to 
meet strength requirements. Repeating the calculations in steps 4.b, 4.c, 
and 4.d for the other span stations and summarizing the results we have the 
results listed in table 3. 

Remarks : 

(1) The maximum number of 0° plies required is four while that for ±45 is 
six. These ply combinations result in an unsymmetric laminate. 

(2) The total number of plies required varies from six in the walls to 
nine in the covers. 

(3) For fabrication convenience, assume constant number of plies through- 
out the box beam. 
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(4) Since the laminate is relatively thin, panel buckling will control 
the design. 

<5) Decide on a laminate configuration and check the panel buckling load. 

(6) A reasonable symmetric laminate configuration is the 10-ply laminate, 
0.050 in. thick, as follows: four plies at 0°, four plies at ±45°, and two 

plies at 90° as was already mentioned. Using conventional notation this 
laminate configuration is designated: 

C±45°/02/90] s 


Step 5: 

Calculate the composite stresses at the four-span stations. These 
stresses are calculated by dividing the in-plane load at that station with the 
laminate thickness (N cxx /tjj_, etc.). Summarizing the results in tabular 
form, we have the results listed in table 4. 

The composite stresses summarized above will be used to check the buckling 
stresses and the ply stresses as described below. 


Step 6: 

Calculate the buckling stresses of the panels at each bay. In order to 
expedite these calculations, we use the midbay panel dimensions and apply the 
stresses summarized above. This approach is reasonable for preliminary design. 
However, it needs to be checked with finite element analysis for more accurate 
results. 


Tabulating the results we have the summary shown in table 5. 


We calculate the buckling stresses by using the approximate interaction 
equation (ref. 1): 


a 

cxx 

(cr) 

a 

CXX 



< 1 


where 


(cr> 

a cxx 


2. 2 C 
* ¥ 


l2b <’ - “cxy'V 



(cr) 

a cxy 


7* 2 t 2 E 


,2b <’ - “cxAyx’ 


(1 < | < 2 ) 


E = 


V 


4E cxx E cyy G cxy 
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The moduli and Poisson's ratios for a [±45/02/90] s AS/E angleplied 
laminate are (ref. 1): 


E C xx = 9.6 mps i ; E C yy = 6.5 mpsi; 

G C xy = 2.3 mpsi ; x y = 0 . 3 3 1 ^cyx = 0*22 

Substituting these values, we calculate 

E = 3 -^4 x 9.6 x 6.5 x 2.3 mpsi = 8.31 mpsi 

First we check the buckling stresses of the top cover at the first midbay 
(0 to 20) 


(cr) 

a 

cxx 


it 2 (0.05) 2 x 8.31 x 10 6 
12(18. 4) 2 (1 - 0.33 x 0.22) 



= 218 psi 


(cr) 

a 

cxy 


— (0 - 05 ^ — x 8 - 31 x = ±380.9 psi 

1 2( 1 8 .4; ( 1 - 0.33 x 0.22) 


Using these buckling stresses in the combined stress interaction equation 
above, we calculate: 


-79 200 ( -6700 V , 

-218 + y-380.9 ) - 

672.7 >> 1 

indicating that the panel will buckle. 

Remark: These panel thicknesses are too low to resist buckling due to the 

applied load stresses in the top cover and back side wall panels since 
all these panels are subjected to combined compressive and shear 
stresses. On the other hand, the panels in the bottom cover and front 
side wall may not buckle because these panels are subjected to com- 
bined tensile and shear stresses. The most critical case is the bot- 
tom cover at the third midbay (40 to 60). The calculated individual 
buckling stresses for this panel are: 

- 7,4 1”'; ‘S’ • 942 P si 

and combined in the interaction equation: 

59 260 ( 27 400 \ 2 . 

-714 + y 942 J - 1 

763 >> 1 and this panel will also buckle. 

The conclusion from the above calculations is that the panel thickness 

sized for strength is too thin to resist buckling. 
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Step 7: 


At this point, we can consider several alternatives to increase the panel 
buckling resistance. The obvious ones are: (1) increase the panel thickness, 
(2) reduce the panel edge dimensions by using inner walls and additional bulk- 
heads, and (3) use combinations of these. 


First we check alternative (1) - increase panel thickness. Calculate 
panel thickness to resist buckling stress. Since the buckling stress varies 
with the thickness squared and assuming panel thickness in multiples of 
[*45/02/90 ] s we calculate a thickness for the compressive stress: 


t 


c 



use 0.45 in. ~ 0.50 in. 


for combined stress. 


This results in a [±45/02/90 ]]qs symmetric laminate with 100 plies. 

This many plies will substantially increase the material and fabrication costs. 

Check buckling stresses for the same panel. 


a cxx = <0.05/0.50) 79 200 psi = 7920 psi 
0 CX y = (0.05/0.50) 6700 psi = 670 psi 

a cxx } = ^-SO/O-OB) 2 218 psi = 21 800 psi 

(rr) ? 

C cxy = (0 - 5/0 - 05) 380.9 P s1 = 38 090 PS i 



Therefore the panel satisfies the combined stress buckling interaction equation 
with a margin of safety 1 - 0.363 = 0.64. 


Alternative (2) will increase the buckling stresses but will not reduce 
the stresses due to applied loads. Alternative (3) on the other hand will 
increase the buckling stresses and also reduce the stresses due to applied 
loads. 


Check buckling stresses by adding an inner vertical wall through the box 
beam center. Since this will reduce the panel edge dimension by 2 it will 
increase the buckling stress by 4. Therefore we assume 60-ply laminate as 
follows: 


9 



[±45/0 2 /90] 6s 


with 0.3 in. thickness. 

The stresses in the panel at the first bay are: 

o cxx = (0.05/0.3) 79 200 psi = 13 200 psi 
o cxy = (0.05/0.3) 6700 psi = 1117 psi 
The corresponding buckling stresses are 

a cxx } = 4(0.3/0.05)2 218 psi = 31 392 psi 

a cxy } = 4(0 - 3/0 -°5) 2 380.9 psi = 54 850 psi 
The combined stress buckling interaction equation is 

13. 200 + (jm V _ 0 42 < 1 OK 

31 392 1 54 850 J ' U * 4 U 

and the margin of safety MOS = 0.58. 

This is a lighter weight design (by 30 percent) compared to that of 
alternative (1). 


Step 8: Summarize Design 

The designed box-beam, therefore is a 60-ply [±45/C>2/90]6s laminate 
0.3 in. thick with an inner vertical wall through the box-beam center. The 
panel geometry at mid bays and the respective stresses at the bulkheads are 
summarized in table VI. 


Step 9: 

(a) Vertical displacement , 
(neglecting box beam weight) 


Check Tip Displacements 
- This displacement is calculated from 


03 = 


P z r 


3E cxx I cyy 
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where Icyy is calculated at midspan as follows: 

,2 


I 


cyy 


“ Ml) *c * 3(fe)h3* c ], 


= ^2 x IsfjS 

= (126.6 + 31 .6)in. 4 


in, 


(0.3) + 3 I 2 ( 7 . 5 ) 3 ( 0 


. 3)Jin.' 


= 158.2 in. 

The values of the variables required in the above equation are 
P z = 13 200 lb 

2. = 60 in. 

E cx* ■ 9 ' 6 "I 551 
'cyy ■ '58-2 in. 4 


•3 


w = 


13 200 1b x 60 x 60 x60 in 
3x9 600 000 (lb/in? ) x 158.2 inf 


w = 0.63 in. < 0.90 (1.5% 60 in.) in. O.K. 
0.90 


MOS = 


0.63 


-1.0= 0.43 


(b) Lateral displacement . - This displacement is calculated from 

,3 


v = 


y 


^cxx'czz 


where, again, I czz is calculated at midspan as follows: 

4 


I 


czz 


(!) 


2(7.5) 


^|^(15) 3 (0.3) i 


= [253.1 + 168.831 n. 
= 421.9 in. 4 


(0.3) + 2(4^|( 1 5) J (0. 3) | in, 
4 
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The values of the variables required to calculate the lateral displacement are 
P y = 6600 lb 

ft = 60 in. 

E cxx ■ 9 - 6 ”P si 
I C2Z - 421.9 in. 4 

= 6600 lb x 60 x 60 x 60 in. 3 

3x9 600 000 (lb/in. 2 ) x 421.9 in. 4 

v = 0.12 in. < 0.90 in. O.K. 

MOS = -1 = 6.5 

(c) Angle of Twist . - This angle is calculated from 

M ft 

9 “ 3S~ 

where 


J *cxx + *cyy 
and 


G = G = G 
cxz cxy 

The values of the variables required to calculate the twist angle are 
M = 200 000 in. -lb 

A 

ft = 60 in. 

J = (158.2 + 421.9) in. 4 = 580.1 in. 4 


G = 2.3 mpsi 
0 = 


200 000 in. -lb x 60 in. 
580.1 in. 4 x 2300 000 lb/in. 2 


rad 


= 0.008994 rad 


9 = 0.52° < 1 .0° O.K. 

M0S ■ 'di - ' - °- 92 
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Step 10: 


Calculate the first flap-wise, edge-wise frequencies, and the first 
torsional frequency. 

(a) The flap-wise (vertical) frequency is calculated from 

1/2 

cyc/sec 



where M is the mass per unit length and I CX x Is the moment of inertia, 
both calculated at midspan (2-covers, 3-walls and 4-bulkheads) 

r 4 1 - 1 

- V'* _- 


M = 


2b + 3h + l £ (bxh) . 

i = l 1 


30 + 22.5 + (478) 


0.3 x 0.06 1U /...2 
386.4 lb/sec 


= 0.0028 


1b sec 1 
in. in, 


The values of the variables required to calculate this frequency are 
9. = 60 in. 

E cxx ■ 9 ' 6 m P s1 
'cxx - ,58 - 2 in - 4 


M = 0.0028 


lb sec 1 


in. in. 


CO_ = 


(1.9 x 1.9) 


z 2ir(60 in. x 60in.) 


9600 000(lb/in. 2 )158.2 in. 4 


0.0028 


lb/sec 1 


in. in. 


1/2 


cyc/sec 


co = 117.5 cyc/sec > 100 cyc/sec O.K. 


MOS = - 1 =0.18 


(b) The edge-wise lateral frequency is calculated from 

“ 1/2 
'll! I 

M 


1 (} .9\ 2 E cxx I cyy 

w y " 2ir V P 


This equation differs from co z only in I C yy. We can expedite the 
calculation. 
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1/2 


Uz) . 9 ^ 

~ \ 1 58 . 2 / 


x 117.5 cyc/sec 


Uy s; 191.9 cyc/sec > 500 cyc/sec O.K. 


M0S = "T^T " 1 = 0,23 


(c) The torsional frequency is calculated from 


u t ~ 42, 


_L fej 


The values for the variables are 
a = 60 in. 

G = 2.3 mpsi 

2 

g = 386.4 in/sec 
p = 0.06 lb/in. 3 


“t ~ 4 x 60 


1 .0 r 2300 000( 1 b/i n . 2 ) x 386.4 in. /sec 2 


0.06 lb/in. 3 


j, 507.1 cyc/sec > 450 cyc/sec O.K. 


M0S = _1 = 0J3 


Step 11: Check Local Panel Vibration 

We calculate the first frequency for the first bay panel (0-20 span) 
assuming a rectangular panel with midside dimensions. This frequency is given 
by 


^UZpU - UcxAyx > 


<’ ♦ ^cyx^cxx * E cyy c4 


+ 4C (1 - u u )G 

cxy cyx cxy 
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where 


t panel thickness (in.) 

a panel x-edge dimension (in.) 

2 

g gravity acceleration (in. /sec ) 

3 

p composite laminate density (Ib/in. ) 

C a/b where b is the panel y-edge dimension (in.) 

where E CX x» Ecyy» anc * Gcxy are the composite laminate moduli and where 
u CX y and \> C y x are composite laminate Poisson's ratios. 

The values for the variables in the frequency calculation are 


= 0.3 in. 

E cxx ■ 9 - 6 n " )Si 

= 20 in. 

E cyy ' 6 ' 5 

2 

= 386.4 in. /sec 

G cxy - 2 3 n, l Jsi 

= 0.06 1 b/i n . 3 

”cxy = °' 33 

= (20/9.2) = 2.17 

v = 0.22 


cyx 


Substituting these values in the frequency equation 


_ ir(0.3) [" 386.4 

“ “ 2 x 20 x 20 12 x 0.06(1 - 0.33 x 0.22) 


1/2 


+ 4(2 . 1 7 ) 2 ( 1 - 0.33 x 0.22) x 2.3 + 6.5 x (2.17) 4 


(1 + 2 x 0.22 x 2.17) x 9.6 
1/2 




x 1000 cyc/sec 


cj = 414.4 cyc/sec 

This frequency is greater than the first two (<d 2 and coy) frequencies of 
the box beam. Therefore, no local vibration will occur prior to the first box 
beam modes. However it could occur prior to the first torsional mode. 


Step 12: Check Ply Stresses 

The ply stresses are determined through the use of the ply stress 
influence coefficients as described in detail in Ref. 1. The ply stress 
influence coefficient for the laminate selected are shown in table VII. 

The ply stress a^n is calculated as follows (include only nonzero 
coefficients) 

°2.1 1 = 1-98 a^xx _ 0.56 a^yy + 9.35 AT + 627 M 
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Examining the ply stress influence coefficients we see that for 

the +45° and 0^2 for the -45° ply have relatively large values. From 
the panel composite stresses summary in Step 8 we check the bottom cover in 
the bays 2 and 3 as follows: (neglecting temperature and moisture). 

(a) +45°-Ply (3-bay bottom cover) 
o£ll = 0.70 o cxx + 4.0 a cxy 

= 0.70 x 9877 + 4.0 x 4567 


= 25182 psi < 


220 000 psi 
220 000 


MOS = 


25182 


O.K. 

-1 = 7.74 


(b) -45°-Ply (1-bay bottom cover) 

CT 2.22 = °- 55 CT cxx + °- 19 a cxy 

= 0.55 x 13200 + 0.19 x 2217 


= 7681 psi < 8000 psi O.K. 


MOS 


8000 _ 
7681 


0.04 


(c) Check the above ply stresses by including residual and moisture 
stresses. For residual stresses AT = -300 °F which is the difference 
between the cure and room temperatures. For the moisture stresses assume 
M = 1 percent by weight. 

+45°-Ply: 

0q,1 1 = 0.70 Qqxx ^ • OO 0 ^ x y + 21.02 AT + 1413 M 

= 0.70 x 9877 + 4.0 x 4567 + 21.02 (-300) + 1413 (1) 

= 6914 + 18268 - 6306 + 1413 


= 20289 psi < 220 000 psi O.K. 


MOS 


220 000 
20289 


9.8 


-45°-Ply: 

a 0.22 = 0-55 a cxx + 0.19 a cxy - 1876 AT - 1263 M 

= 0.55 x 13200 + 0.19 x 2217 - 18.76 (-300) - 1263 (1) 
= 7620 + 421 + 5628 - 1263 
= 12046 psi > 8000 psi N.G. 
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MOS 


8000 

12046 


-0.34 


This last calculation indicates that the -45°-Ply will crack in transverse 
tension due to combined design mechanical and environmental loads. The last 
calculation also illustrates the significance of residual stresses in compos- 
ites. Since this is a matrix failure mode, we check the ply stress' at speci- 
fied mechanical loads. Recall that the design loads are two times the 
specified loads. 



Step 13: Summarize Design Results 

(a) Laminate configuration [±45/02/90]gs 0.3 in. thick 

(b) Box-beam design - uniform laminate thickness, two intermediate 
bulkheads, and one inner wall located at the box beam center (see fig. 2) 

(c) Box-beam weight = 66 lb (composite volume times density) 

(d) Tip displacements MOS 

Displacement MOS 
w 0.43 

v 6.50 


0 0.92 

(e) Buckling load MOS = 0.58 

(f) Vibration frequencies MOS 


Box beam 

MOS 

frequency 


u z 

0.18 

COy 

0.23 

co t 

0.13 
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Panel 1-bay 
top cover 1.16 


(g) Ply stresses 

+45°-P 1 y (3-bay-bottom cover) 


MOS 

longitudinal stress 

mechanical loads 7.7 

mechanical and environmental 9.8 

-45°-P 1 y (1-bay-bottom cover) 

transverse stress 

mechanical loads 0.04 

mechanical and environmental -0.34 


CONCLUDING REMARKS 

Step-by-step design procedures are described which can be used for the 
preliminary design of composite box beams. The various calculations in these 
procedures are arranged so that they can be performed using a pocket calcula- 
tor. The sample calculations are for the design of a cantilevered composite 
box beam subjected to end loads. The composite laminate is selected to satisfy 
design requirements for local buckling, tip displacements, beam and panel 
vibrations, and ply stresses including thermal and hygral (moisture) stresses. 
The procedures and the sample calculations illustrated can be used for the 
preliminary design of composite built-up structures in general. 
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TABLE I. - TYPICAL PROPERTIES OF UNIDIRECTIONAL COMPOSITES AT ROOM TEMPERATURE 


Properties 

Symbol 

Units 

Boron/ 

epoxy 

Boron/ 

poly 

mide 

S-glass 

/epoxy 

Modmor 1/ 
epoxy 

Modmor 1/ 
polymide 

Thornel 

300/ 

epoxy 

Kevlar 

49/ 

epoxy 

Graphite 

AS/epoxy 

Fiber volume ratio 

kf 


0.50 

0.49 

0.72 

0.45 

0.45 

0.70 

0.54 

0.60 

Density 

n 

lb/in 3 

0.073 

0.072 

0.077 

0.056 

0.056 

0.058 

0.049 

0.057 

Longitudinal thermal 

nil 

10" 6 In 

3.4 

2.7 

2.1 


0.0 

0.01 

-1.60 

0.40 

coefficient 


/in/ *F 









Transverse thermal 

n22 

10- 6 in 

16.9 

15.8 

9.3 

18.5 

14.1 

12.5 

31.3 

16.4 

coefficient 


/in/ *F 









Longitudinal modulus 

E iii 

10 6 psi 

29.2 

32.1 

8.8 

27.5 

31.3 

21.0 

12.2 

16.0 

Transverse modulus 

Et22 

10 6 psi 

3.15 

2.1 

3.6 

1.03 

0.72 

1.5 

0.70 

2.2 

Shear modulus 

G ll2 

10 6 psi 

0.78 

1.11 

1.74 

0.9 

0.65 

1.0 

0.41 

0.72 

Major Poissons's ratio 

m 2 


0.17 

0.16 

0.23 

0.10 

0.25 

0.28 

0.32 

0.25 

Minor Poissons's ratio 

v l21 


0.02 

0.02 

0.09 


0.02 

0.01 

0.02 

0.34 

Longitudinal tensile 

s illT 

psi 

199 000 

151 000 

187 000 

122 000 

117 000 

218 000 

172 000 

220 000 

strength 











Longitudinal compres- 

s tllC 

psi 

232 000 i 

158 000 

119 000 

128 000 

94 500 

247 000 

42 000 

180 000 

sive strength 











Transverse tensile 

s i22T 

psi 

8100 

1600 

6670 

6070 

Z150 

5850 

1600 

8000 

strength 











Transverse compres- 

s i22C 

psi 

17 900 

9100 

23 500 

28 500 

10 200 

35 700 

9400 

36 000 

sive strength 











Intralaminar shear 

$tl2S 

psi 

9100 

3750 

6500 

8900 

3150 

9800 

4000 

10 000 

strength 











Longitudinal moisture 

6f.ll 

1CT 2 in 

0.003 

0.003 

0.014 

0.003 

0.003 

0.006 

0.008 

0.006 

coefficient 











Transverse moisture 

B*22 

ID" 2 in 

0.168 

0.168 

0.128 

0.129 

0.129 

0.129 

0.151 

0.129 

coefficient 











Glass transition 

t GD 

*F 

420 

700 

420 

420 

700 

420 

420 

420 

temperature 











(estimate) 












TABLE II. - FORCES IN COVERS AND SIDEWALLS 


Span 

station, 
i n . 

Length, 
i n . 

Covers 

Side Walls 

Depth, 

in. 

Width, 
i n . 

N cxx ’ 

ki s 

N 

lN cxy » 

ki s 

Oepth, 
i n . 

D 

N 

cyy 1 

ki s 

M CXZ 1 

ki s 

top 

bot 

top 

bot 

front 

back 

front 

back 

■m 

60 

10 

20 


3960 

-335 

665 

20 

10 


-1980 

160 

1160 

■ 

40 

8.3 

16.7 

-3809 

3809 

-523 

920 

16.7 

8.3 


-1904 

74 

1516 

HiKiifl 

20 

6.7 

13.3 

-2963 

2963 

-874 

1370 

13.3 

6.7 


-1481 

-137 

2107 

4. 60 

0 

5 

10 


— 

-2330 

+2330 

10 

5 



-3320 

3320 
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TABLE III. - PLIES IN COVERS AND SIDEWALLS 


Span 

station 

Plies in the covers 

Plies in 

the wal 1 s 

0° 

±45° 

90° 

Total 

0° 

+45° 

90° 

Total 

1 . 0 

4 

mm 

2 

8 

2 


2 

6 

2. 20 

4 

mm 

2 

8 

2 


2 

6 

3. 40 

3 

mm 

2 


2 


2 

8 

4. 60 

- 

H 

- 


- 

6 

- 

6 


TABLE IV. - PLY STRESSES IN COVERS AND WALLS 


Span 
stati on 

1 . 0 
2. 20 

3. 40 

4. 60 

Covers (psi ) 

Side walls (psi ) 

Top 

Bottom 

Front 

Back 

a cxx 

a cxy 

a cxx 

a cxy 

°cxx 

a cxz 

a cxx 

CT cxz 

-79 200 
-76 180 
-59 260 

- 6 700 
-10 460 
-17 480 
46 600 

79 200 
76 180 
59 260 

13 300 
18 400 
27 400 
46 600 

39 600 
38 080 
29 600 

3 200 
1 480 
- 2 740 
-66 400 

-39 600 
-38 080 
-29 600 

23 200 
30 320 
42 140 
66 400 


TABLE V. - BUCKLING STRESSES 


wamm 

m 

Bay/span station 


ii 

(■I 

B 

1 (0-20) 

2 (20-40) 

3 (40-60) 










■ 

M 

■ 

Covers Walls 

Covers Walls 

Covers Walls 

Geometry, 

a 

in. 

20.0 20.0 

20.0 20.0 

20.0 20.0 


3 


18.4 9.1 

15.0 7.5 

11.6 5.8 


^c 


.050 .050 

.050 .050 

.050 .050 

Stresses, psi 
top cover 

-79 200 

-76 180 

-59 260 


7 cxy 


-6 700 

-10 460 

-17 480 

Bottom cover 

79 200 

76 180 

59 260 


7 cxy 


13 300 

18 400 

27 400 

Side walls 
front 

^cxx 


39 600 

i 38 080 

29 600 

c 

T cxz 


3 200 

14 801 

- 2 740 

back 

a cxx 


-39 600 

-38 080 

-29 620 

a cxz 


23 200 

| 30 320 

42 140 
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TABLE VI. - FINAL DESIGN STRESSES 


Bay/span station 


1 ( 0 - 20 ) 


2 (20-40) 


3 (40-60) 


Geometry, 

in. 


Stresses, psi 
Top cover 


Bottom cover o r 


Walls 

front 


zers Walls 

Covers Walls 

Covers Walls 

20 20 

20 20 

20 20 

9.2 4.6 

7.5 3.8 

5.8 2.9 

0.3 0.3 

0.3 0.3 

0.3 0.3 

-13200 

-12697 

-9877 

- 1167 

-17433 

-2913 

13200 

12697 

9877 

2217 

3067 

4567 

6600 

6347 

4933 

-200 

- 938 

-2274 

-6600 

-6347 

-4933 

3133 

1161 

5240 

1467 

1767 

2189 



TABLE VII. - COMBINED HYGROTHERMOMECHANICAL LOAD STRESS ASSESSMEMT 


Ply/ply 

stress 




Ply Stress Influence 

Coefficients 

Per unit 


Per unit 
composi te 

Per unit 

composite moisture, 

compo 

a cxx 

site st 
°cyy 

ress 

a cxy 

temperature, 

(°F) 

( 1 percent by wei ght) 


-0.56 

0 

9.35 

627 


.18 

0 

-18.26 

-1227 

0 

0 

0.27 

0 

0 

0.70 

1.25 

4.00 

21.02 

1413 

0.55 

0.10 

-0.19 

-18.76 

-1263 

-0.09 

0.12 

0 

0.75 

53 

0.70 

1.25 

-4.00 

21.02 

1413 

0.55 

0.10 

0.19 

-18.76 

-1263 

0.09 

-0.12 

0 

-.75 

-53 

-0.57 

3.05 

0 

32.69 

2198 

0.12 

0.01 

0 

-19.37 

-1298 

0 

0 

-0.27 

0 

0 
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